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Abstract: We introduce a new jet shape — A^-subjettiness — designed to identify boosted 
hadronically-decaying objects like electroweak bosons and top quarks. Combined with a 
jet invariant mass cut, iV-subjettiness is an effective discriminating variable for tagging 
boosted objects and rejecting the background of QCD jets with large invariant mass. In 
efficiency studies of boosted W bosons and top quarks, we find tagging efficiencies of 30% 
are achievable with fake rates of 1%. We also consider the discovery potential for new 
heavy resonances that decay to pairs of boosted objects, and find significant improvements 
are possible using A^-subjettiness. In this way, iV-subjettiness combines the advantages of 
jet shapes with the discriminating power seen in previous jet substructure algorithms. 
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1. Introduction 

The Large Hadron Collider (LHC) will search for new physics by probing a previously 
unexplored kinematic regime. Most new physics scenarios that provide a solution to the 
hierarchy problem predict that the LHC will produce new heavy particles with decay chan- 
nels involving top quarks, W / Z bosons, and Higgs bosons. In addition, many extensions 
of the standard model including technicolor and Higgs compositeness invoke new heavy 
resonances within the LHC reach with large branching fractions to pairs of gauge bosons 
and top quarks. Therefore, a key task in the search for physics beyond the standard model 
is to efficiently identify final state electroweak gauge bosons and top quarks in a variety of 
kinematic configurations. 

With its current center-of-mass energy of y/s = 7 TeV, the LHC is already able to 
produce new TeV-scale resonances which can decay to highly boosted electroweak bosons 
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and/or top quarks. For a large enough boost factor, the decay and fragmentation of such a 
boosted object yields a collimated spray of hadrons which a standard jet algorithm would 
reconstruct as a single jet. Thus, standard reconstruction methods for electroweak bosons 
and top quarks become ineffective due to the immensely large background of ordinary QCD 
jets. One possibility is to focus on channels where the boosted object decays leptonically, 
though such methods discard much of the original signal and may therefore not be optimal 
for detecting new heavy resonances. 

Recently, there has been considerable progress in identifying boosted hadronically- 
decaying objects using jet substructure techniques. Algorithmic methods use information 
from the jet clustering procedure to extract the internal structure of jets ^ ||, [|, || |(| , 
and are able to successfully distinguish between jets originating from boosted electroweak 
boson and top quarks (denoted here as U W jets", "top jets", etc.) and those originating 
from light quarks or gluons ("QCD jets"). Jet shape methods efficiently tag boosted 
objects with jet-based observables that take advantage of the different energy flow in the 
decay pattern of signal jets and background jets. In addition, there are "jet grooming" 
techniques such as filtering || pruning [10, pjj] , trimming jE^], and their combinations 
Jl3| ] which aid in the identification of boosted objects by reducing the smearing effects of 
jet contamination from initial state radiation, underlying event activity, and pileup. Taken 
together, these jet substructure methods show much promise for enhancing searches for 
new physics in all-hadronic decay channels. 1 

In this paper, we introduce a new tagging method for boosted objects based on a 
novel jet shape dubbed "iV-subjettiness" and denoted by t^. This variable is adapted 
from the event shape "iV-jettiness" advocated in Ref. [28] as a way to veto additional jet 
emissions and define an exclusive jet cross section. Here, we take advantage of the multi- 
body kinematics in the decay pattern of boosted hadronic objects, and use iV-subjettiness 
to effectively "count" the number of sub jets in a given jet. We find that ti/t\ is an effective 
discriminating variable to identify two-prong objects like boosted W, Z, and Higgs bosons, 
and T3/V2 is effective for three-prong objects like boosted top quarks. 

Compared to previous jet substructure techniques, iV-subjettiness has a number of ad- 
vantages. First, to identify boosted objects, one wants to find jets that contain two or more 
lobes of energy. While previous jet shape measures do capture the deviation of a jet from 
a one-lobe QCD-like configuration, TV-sub jettiness is a more direct measure of the desired 
energy flow properties. Second, it is convenient to be able to adjust the relative degree of 
signal efficiency and background rejection without having to perform computationally in- 
tensive algorithmic adjustments. Like for other jet shape methods, tjv can be calculated for 
every jet, and a flexible one-dimensional cut on a function f(ri, . . . , r^r) can determine the 
efficiency /rejection curve. Similarly, the set of tn values can be used as input to a multi- 
variate discriminant method for further optimization. Third, iV-subjettiness is an inclusive 
jet shape and can be defined and calculated without reference to a jet substructure algo- 
rithm. 2 This will likely make iV-subjettiness more amenable to higher-order perturbative 



Additional recent related work appears in Refs. ju], |fj, ^ |TJ [TJ, [TJ, |(J [H], f2§ |2§ ||, |§ & 
2 As we will discuss below, for computational purposes, we use a definition of iV-subjettiness that still 
has residual dependence on a clustering procedure. See Sec. U/2 for further discussion. 
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calculations and resummation techniques (see, e.g. recent work in Ref. |2^, |3(|) compared 
to algorithmic methods for studying substructure. Finally, iV-subjettiness gives favorable 
efficiency/rejection curves compared to other jet substructure methods. While a detailed 
comparison to other methods is beyond the scope of this work, we are encouraged by these 
preliminary results. 

The remainder of this paper is organized as follows. In Sec. |2|, we define ./V-subjettiness 
and discuss some of its properties. We present tagging efficiency studies in Sec. ||, where we 
use iV-subjettiness to identify individual hadronic W bosons and top quarks, and compare 
our method against the YSplitter technique [§, || U and the Johns Hopkins Top Tagger || . 
We then apply iV-subjettiness in Sec. || to reconstruct hypothetical heavy resonances de- 
caying to pairs of boosted objects. Our conclusions follow in Sec. ||, and further information 
appears in the appendices. 

2. Boosted Objects and iV-subjettiness 

Boosted hadronic objects have a fundamentally different energy pattern than QCD jets 
of comparable invariant mass. For concreteness, we will consider the case of a boosted 
W boson as shown in Fig. [l], though a similar discussion holds for boosted top quarks or 
new physics objects. Since the W decays to two quarks, a single jet containing a boosted 
W boson should be composed of two distinct — but not necessarily easily resolved — hard 
subjets with a combined invariant mass of around 80 GeV. A boosted QCD jet with an 
invariant mass of 80 GeV usually originates from a single hard parton and acquires mass 
through large angle soft splittings. We want to exploit this difference in expected energy 
flow to differentiate between these two types of jets by "counting" the number of hard lobes 
of energy within a jet. 

2.1 Introducing iV-subjettiness 

We start by defining an inclusive jet shape called "./V-subjettiness" and denoted by r^r. 
First, one reconstructs a candidate W jet using some jet algorithm. Then, one identifies 
N candidate subjets using a procedure to be specified in Sec. |2.2| . With these candidate 
subjets in hand, tn is calculated via 

tn = -r y^PT,fcmin{Ai?i, fc , A-^.fc, • • • ' AR N,k}- (2-1) 
d °k 

Here, h runs over the constituent particles in a given jet, pr,k are their transverse momenta, 
and ARj^ = \J [At]) 2 + (A0) 2 is the distance in the rapidity-azimuth plane between a 
candidate subjet J and a constituent particle k. The normalization factor do is taken as 

do = ^2pT,kRo, (2.2) 

k 

where Ro is the characteristic jet radius used in the original jet clustering algorithm. 

It is straightforward to see that tn quantifies how iV-subjetty a particular jet is, or 
in other words, to what degree it can be regarded as a jet composed of N subjets. Jets 
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Boosted W Jet, R = 0.6 




Figure 1: Left: Schematic of the fully hadronic decay sequences in (a) W + W~ and (c) dijet QCD 
events. Whereas a W jet is typically composed of two distinct lobes of energy, a QCD jet acquires 
invariant mass through multiple splittings. Right: Typical event displays for (b) W jets and (d) 
QCD jets with invariant mass near mw The jets are clustered with the anti-fey jet algorithm ]3l| ] 
using R = 0.6, with the dashed line giving the approximate boundary of the jet. The marker size 
for each calorimeter cell is proportional to the logarithm of the particle energies in the cell. The 
cells are colored according to how the exclusive kx algorithm divides the cells into two candidate 
subjets. The open square indicates the total jet direction and the open circles indicate the two 
subjet directions. The discriminating variable t%/tx measures the relative alignment of the jet 
energy along the open circles compared to the open square. 



with ttv ~ have all their radiation aligned with the candidate subjet directions and 
therefore have ./V (or fewer) subjets. Jets with t^t 3> have a large fraction of their energy 
distributed away from the candidate subjet directions and therefore have at least N + 1 
subjets. Plots of t\ and T2 comparing W jets and QCD jets are shown in Fig. |2|. 

Less obvious is how best to use tat for identifying boosted W bosons. While one might 
naively expect that an event with small T2 would be more likely to be a W jet, observe that 
QCD jet can also have small T2, as shown in Fig. 2(b) . Similarly, though W jets are likely 
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65 GeV < m. < 95 GeV 65 GeV < m, < 95 GeV 

i I 




(a) (b) 

Figure 2: Distributions of (a) t\ and (b) T2 for boosted W and QCD jets. For these plots, we 
impose an invariant mass window of 65 GeV < mj Ct < 95 GeV on jets of R = 0.6, pr > 300 GeV, 
and \t}\ < 1.3. By themselves, the tn do not offer that much discriminating power for boosted 
objects beyond the invariant mass cut. 
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(a) (b) 

Figure 3: (a): Distribution of T2/T1 for boosted W and QCD jets. The selection criteria are the 
same as in Fig. |[ One sees that the T2/T1 ratio gives considerable separation between W jets and 
QCD jets beyond the invariant mass cut. (b): Density plot in the T1-T2 plane. Marker sizes are 
proportional to the number of jets in a given bin. In principle, a multivariate cut in the t\-ti plane 
would give further distinguishing power. 



to have large t±, QCD jets with a diffuse spray of large angle radiation can also have large 
Ti, as shown in Fig. 2(a). However, those QCD jets with large t\ typically have large values 



of T2 as well, so it is in fact the ratio T2/T1 which is the preferred discriminating variable. 
As seen in Fig. [3(a) , W jets have smaller T2/T1 values than QCD jets. Of course, one can 
also use the full set of tjv values in a multivariate analysis, as suggested by Fig. 3(b) , and 
we will briefly explore this possibility in Sec. 3.4. 

As mentioned in the introduction, iV-subjettiness is adapted from the similar quantity 
iV-jettiness introduced in Ref. |28|]. There are three important differences: the sum over 
k only runs over the hadrons in a particular jet and not over the entire event, we do not 
have candidate (sub)jets corresponding to the beam directions, and our distance measure 
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is only longitudinally boost invariant and not fully Lorentz invariant. The definition of 
tn is by no means unique, and some variations are discussed in App. [A], though Eq. ( |2.1| ) 
appears to be well-suited for boosted object identification. 

2.2 Finding Candidate Subjets 

A key step for defining A-subjettiness is to appropriately choose the candidate subjet 
directions. As also mentioned in Ref. |28|], ideally one would determine rjv by minimizing 
over all possible candidate subjet directions, analogously to how the event shape thrust is 
defined |32|. In that case, tn is a strictly decreasing function of A, and < tn/tn-i < 1. 

In practice, such a minimization step is computationally intensive, and for the studies 
presented in this paper, we will determine the candidate subjet directions by using the 



exclusive- kjp clustering algorithm 33, 34], forcing it to return exactly A jets. This algo- 
rithmic method for finding the candidate subjet momenta is not perfect, and though not 
displayed in Fig. 3(a), a small fraction («C 1%) of jets have tn > tjv_i, a "feature" which 



would be eliminated with a proper minimization procedure. That said, we have not found 
many cases where the exclusive Jzt jet algorithm identifies completely incorrect subjets, 
and the additional event displays in App. |B] show that the exclusive kx algorithm finds 
subjets that look reasonable. 

Once the candidate subjets are identified, A-subjettiness is a proper inclusive jet shape. 
Since Eq. ( |2.lD is linear in each of the constituent particle momenta, tn is an infrared- 
and collinear-safe observable. That is, the addition of infinitesimally soft particles does 
not change A-subjettiness (infrared safety), and the linear dependence on the particle 
momenta combined with the smooth angular dependence ensures that the same tm is 
obtained for collinear splittings (collinear safety). Crucially, the candidate subjets used in 
A-subjettiness must be determined via a method that is also infrared- and collinear-safe, 
something that is automatic with a minimization procedure or by using declustering. 

2.3 Summary 

To summarize, .A-subjettiness is an inclusive jet shape that offers a direct measure of 
how well jet energy is aligned into subjets, and is therefore an excellent starting point for 
boosted object identification. The ratio tn/tn_\ is an easily adjustable offline parameter 
which can be varied to adjust signal efficiency/background rejection without having to redo 
the clustering of the particles in an event. While there is some residual jet algorithm depen- 
dence in the identification of the original seed jet and in the identification of the candidate 
subjets, this latter effect could be completely removed by using a minimization procedure 
at the expense of introducing more computational complexity. Though we will not attempt 
to do so here, we suspect that A-subjettiness will lend itself better to theoretical studies 
than algorithmic boosted object tagging methods, either in fixed-order or resummed QCD 
calculations. As we will see in the next two sections, .A-subjettiness compares favorably to 
other boosted object tagging methods in terms of discriminating power. 

Finally, in the above discussion, we used boosted W bosons just as an example, and 
a similar discussion holds for Z bosons and Higgs bosons. Also, A-subjettiness will be 
effective for identifying boosted top quarks. A top quark with mass of 175 GeV decays to 
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Boosted Top Jet, R = 0.8 




Figure 4: Left: Decay sequences in (a) tt and (c) dijet QCD events. Right: Event displays for 
(b) top jets and (d) QCD jets with invariant mass near mtop- The labeling is similar to Fig. ^, 
though here we take R = 0.8, and the cells are colored according to how the jet is divided into 
three candidate subjets. The open square indicates the total jet direction, the open circles indicate 
the two sub jet directions, and the crosses indicate the three sub jet directions. The discriminating 
variable t^Jti measures the relative alignment of the jet energy along the crosses compared to the 
open circles. 

a b jet and a W boson, and if the W boson decays hadronically into two quarks, the top jet 
will have three lobes of energy. Thus, instead of one expects T3/T2 to be an effective 

discriminating variable for top jets. This is indeed the sketched in Figs. f|, [|, |(| 

and @. 
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Figure 5: Distributions of (a) t%, (b) t 2 and (c) T3 for boosted top and QCD jets. For these plots, 
we impose an invariant mass window of 145 GeV < m- ]ct < 205 GeV on jets with R = 0.8, pr > 300 
GeV and \t}\ < 1.3. 
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Figure 6: Distributions of (a) T2/71 and (b) T3/T2 for boosted top and QCD jets. The selection 
criteria are the same as in Fig. ||. We see that T3/T2 is a good discriminating variable between 
top jets and QCD jets. In this paper, we do not explore T2/T1 for top jets, though it does contain 
additional information. 
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Figure 7: Density plots in the (a) ti~t 2 plane and (b) t 2 -t 3 plane for boosted top and QCD jets. 
The selection criteria are the same as in Fig. |^. These plots suggest further improvement in boosted 
top identification is possible with a multivariate method. 
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3. Efficiency Studies 



In this section, we investigate the tagging efficiencies for individual W jets and top jets 
and the mistagging rates for QCD jets. The W boson case is representative of hadronically 
decaying Z bosons and Higgs bosons, though for simplicity we only show W jet results 
in this paper. Here, we select candidate boosted objects using an invariant mass cut 
augmented with an iV-subjettiness criterium and compare our results to the YSplitter 
method @, ||, |j] as well as to the Johns Hopkins Top Tagger §. This study will lay the 
groundwork for the case study in Sec. |] for reconstructing a new physics resonance which 
decay to pairs of gauge bosons or top quarks. 

3.1 Analysis Overview 

The basic criteria for tagging a boosted W boson or top quark is that the jet invariant 
mass should fall near mw or nttop, respectively. For concreteness, we consider the mass 
windows of 65 GeV < m jct < 95 GeV for W jets and 145 GeV < m jet < 205 GeV for top 
jets. We then apply a cut on the r^jrx ratio (for W jets) or the T3/V2 ratio (for top jets), 
where the cut is adjusted to change the relative signal tagging efficiency and background 
mistagging rate. 

For all of our studies, we generate events for pp collisions at a center-of-mass energy 
of yfs = 7 TeV with Pythia 8.135 and perform jet clustering with Fast Jet 

2.4.2 [37, 38]. For the QCD background of light quarks and gluons, we use the default 
QCD dijet production routines. For the boosted W boson signal we use standard model 
W + W~ production, and for the boosted top signal we use standard model ti production. 
In both signal samples, we force the W bosons to decay entirely hadronically, eliminating 
the leptonic decays of the W bosons for which our method is not applicable. Everywhere 
we include multiple interactions, initial- and final-state radiation (ISR/FSR), and hadron 
level decay, though not the effects of event pileup. 

We apply a global analysis cut to isolate boosted central jets, namely jets with trans- 
verse momentum pt > 300 GeV and rapidity \rj\ < 1.3), and we take only the hardest 
jet in each event. 3 In some of our subsequent analyses, we apply harder cuts beyond the 
Pt > 300 GeV restriction, because for more moderately boosted jets, our method becomes 
less effective. 

To partially simulate detector effects and to speed up jet reconstruction, observable 
final-state particles with \r]\ < 4 are collected into "calorimeter" cells arranged in a rectan- 
gular lattice with 80 rapidity (77) and 64 azimuth ((f)) bins (corresponding to approximately 
0.1 x 0.1 sized cells). The calorimeter momenta are interpreted as massless pseudo-particles 
with energy given by the calorimeter energy. 



The seed jets for analysis are determined using the anti-/cr jet algorithm [31|, with 



various jet radii R. To compute r^v, the seed jets are reclustered with the exclusive- £>r 



To reduce computing time for the simulations, we utilize a parton-level momentum phase space cut 
of pr > 200 GeV on all partons studied in this paper. This cut applies both to gluons and light quarks, 
as well as W bosons and top quarks. Cross-checks show that outside this kinematic regime, a negligible 
number of events contain jets with pr > 300 GeV, so no significant bias is introduced. 
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algorithm 34 1 into exactly N candidate subjets, and these subjets are used as input 
to Eq. fl2.ip . Note that Eq. (|2.1[ ) only depends on the three-momenta of the candidate 
subjets. 

Finally, in our study, we used the default Pythia showering algorithm to describe 
ISR/FSR and did not include possible matrix element corrections. While a proper modeling 
of perturbative radiation would certainly affect the specific values of A-subjettiness as well 
as the jet invariant mass distribution, we found that the ratio tn/tn-i is reasonably robust 
to changes in the shower model. Thus, we suspect that using matched multi-jet samples 



1 39(1 would not significantly affect the tagging efficiencies of our A-subjettiness cuts, though 



a detailed study is beyond the scope of this work. 
3.2 Comparing to Other Methods 

To evaluate the performance of A-subjettiness, we want to compare it to previous jet 
substructure methods. The most natural comparison is to the YSplitter technique || |4|. 
In YSplitter, a jet is declustered using the exclusive-^ jet algorithm, and the un,n+i 
scale is the square-root of the kx distance measure at which the jet declusters from A 
subjets to A + 1 subjets. Especially since we are using the exclusive- £>r jet algorithm in 
this paper to define ttv, one might naively think that tn/tn-i and yN,N+\/VN-i,N should 
have similar discriminating power. However, we will see that tn does provide additional 
information beyond the Un,n+i variable. This makes sense, as Vn,n+i measures the "scale" 
at which subjets decompose, whereas tn /tjv-i measures the degree to which jet radiation 
is collimated around the candidate subjet directions. In other words, tn is a more direct 
measure of how A-subjetty a jet really is. 

While a detailed comparison of A-subjettiness and YSplitter is beyond the scope of 
this work, we find that A-subjettiness does compare favorably to YSplitter in a simple, 
non-optimized test. For this naive comparison with the YSplitter method, we use the same 
anti-fey jet clustering parameters and the same invariant mass windows for W jets and 
top jets. We then place cuts on the ratio 2/23/2/12 (for W jets) and 2/34/2/23 (for top jets). 
This is certainly not an ideal use of the YSplitter variables, and we consider potential 
optimizations in Sec. 3A. Plots of various ]Jn,n+i distributions appear in App. |C]. 



Finally, for top jets, we also do one benchmark comparison to the Johns Hopkins Top 
Tagging (JHTT) method, using settings similar to those in Ref. ||. Here, jet clustering 
is performed with the Cambridge- Aachen algorithm [40, 41] combined with the Fast Jet 



JHTT plugin [37] with a jet radius of R = 0.8, 5 P = 0.10, and 5 r = 0.19. As above, we 
only consider the hardest jet in each event, provided it has px > 300 GeV and \r]\ < 1.3. 
We impose the top mass window (145 GeV < mj 0t < 205 GeV) on the entire jet and also 
demand that any two subjets in the JHTT declustering sequence have an invariant mass 
near mw (65 GeV to 95 GeV). Finally, we selected on the W helicity angle by requiring 
cos9h < 0.7. Again, this is not an optimal comparison, but serves as a useful point of 
reference to gauge the effectiveness of A-subjettiness. 
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W jets, R = 0.6, p T > 300 GeV, \r\\ < 1.3 
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Figure 8: Invariant mass distributions for (a) W jets and (b) QCD jets as the iV-subjettiness cut 
on t 2 /ti is varied. Here, the jet radius is R — 0.6 and the jets satisfy px > 300 GeV and < 1.3. 
Since the QCD jet rate decreases faster than the W jet rate as the t 2 /ti cut is tightened, t 2 /ti is 
an effective discriminating variable. 
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.51 


.085 


r 2 /n < 0.3 
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.024 


t 2 /ti < 0.2 


.16 


.0044 


.16 
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.16 : 


.0052 


.15 


.0036 


.16 


.0034 


1% mistag rate 


.26 


.010 


.24 


.010 


.23 : 


.010 


.24 


.010 


.26 


.010 


40% tag efficiency 


.40 


.025 


.40 


.025 


.40 : 


.028 


.40 


.036 


.40 


.045 



Table 1: Tagging efficiencies vs. mistagging rates for W jets : QCD jets with R = 0.6. The top 
row corresponds to just applying the myy invariant mass window (65 GeV to 95 GeV) criteria, and 
the subsequent rows include an additional T2/T1 cut. The bottom two rows indicate the tagging 
efficiencies achievable with a fixed mistagging rate of 1%, and the mistagging rate achievable with 
a fixed tagging efficiency of 40%. 



3.3 Boosted W and Top Results 

We are now ready to investigate the efficiency of the iV-subjettiness method applied to 
W jets and top jets. We start with W jets. In Fig. ^, the invariant mass distribu- 
tions for boosted W jets and QCD jets are shown for boosted central jets with R = 0.6, 
Pt > 300 GeV and \rj\ < 1.3. As the cut on tz/tx is tightened, the boosted W jet sig- 
nal experiences a modest decrease, but the peak of the distribution stays well-centered 
on the W mass. The QCD background decreases by a much greater degree, with no sig- 
nificant sculpting of the distribution. Thus we confirm that t^Jtx is effective as a W jet 
discriminating variable. 

In Table [l], we show the quantitative effect of various t^/tx cuts for different jet pr bins. 
First, note that the mw invariant mass window (65 GeV < m- ]e t < 95 GeV) already acts 
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Figure 9: W jet signal efficiency/background rejection plots. Here, the t<z/t\ cut is varied for 
different (a) jet transverse momenta and (b) jet radii. Unless otherwise indicated, the jets have 
R = 0.6 and 450 GeV < px < 600 GeV. The rightmost points in each plot are for the raw invariant 
mass window criterium alone, and points to the left of these are obtained with additional cuts on the 
T2/T1 ratio. Figure (c) shows a naive comparison against the YSplitter method, where the purple 
lines indicate varying cuts on t^/tx and 2/23/2/12, respectively, while the shaded bands indicate the 



modified linear cuts of Sec. 3.4. In this non-optimized test, ./V-subjettiness compares favorably to 
YSplitter. 




Figure 10: Improvement compared to an invariant mass cut alone for W jets in (a) signal over 
background (S/B) and (b) signal over square- root background (S/V13) by using the t^/tx cut. The 
value of unity on the right corresponds to using only the niw invariant mass window. 



as a useful discriminating variable, with a tagging efficiency of around 70% for mistagging 
rates between 14% and 25% across the whole kinematic range. (The higher mistagging rate 
for more boosted jets occurs because the invariant mass of QCD jets increases with their 
transverse momenta.) By including a modest cut on T2/T1, we can maintain nearly the 
same signal efficiency while halving the mistagging rate. By tightening the T2/T1 cut, we 
can control the degree of signal efficiency and background rejection, and across the whole 
Pt range, a tagging efficiency of 25% can be achieved for a mistagging rate of only 1%. 

In Fig. [|, we plot tagging efficiency curves. These show the effects of the Tijr\ cut in 
different kinematical regimes and for different jet radii. In Fig. 9(b) , the signal efficiency is 
lower for R = 1.0 as compared to R = 0.8, because a larger fraction of W jets fall outside 
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Top jets, R = 0.8, > 500 GeV, |q| < 1 .3 QCD jets, R = 0.8, p T > 500 GeV, |r|| < 1 .3 




V T 2 <0 - 6 



50 100 150 200 250 50 100 150 200 250 

Invariant mass (GeV) Invariant mass (GeV) 

(a) (b) 

Figure 11: Invariant mass distributions for (a) top jets and (b) QCD jets as the iV-subjettiness cut 
on t 3 /t 2 is varied. Here, the jet radius is R = 0.8 and the jets satisfy pr > 500 GeV and \r]\ < 1.3. 
Since the QCD jet rate decreases faster than the top jet rate as the T3/T2 cut is tightened, T3/T2 is 
an effective discriminating variable. Note also that the secondary peak at mw decreases when the 
t 3 /t 2 cut is applied. 



of the mass window due to fact that more ISR is captured in jets with larger radii. As 
a result, the signal efficiency using only a mass cut is significantly reduced. Beyond this 
mass cut, however, our method is relatively insensitive to changes in jet radius, as the slope 
of the efficiency curves does not change considerably with changes in R. In Fig. B7c], we 



compare a T2/V1 cut to a 2/23/2/12 cut, which shows that an iV-subjettiness cut compares 
favorably to a naive YSplitter cut. In Fig. [l(| we show how the signal to background ratio 
(S/B) and signal to square-root background ratio (S/VB) improve as the T2/V1 cut is 
tightened. Compared to only having an invariant mass cut, the S/B ratio can improve by 
as much as an order of magnitude with an improvement of S/^/B of around 50%. These 
improvements will be relevant for the resonance study in Sec. ^j. 

We now consider the analogous tables and plots for top jets, where we use a jet radius 
of R = 0.8, the m top invariant mass window (145 GeV < mj et < 205 GeV), and the 
discriminating variable T3/T2. In Fig. 11, we show the invariant mass distributions for top 
jets and QCD jets, considering px > 500 GeV for both. Note that the T3/T2 cut decreases 
the background rate faster than the signal rate. In addition to the top mass peak, there 
is a secondary peak at m\y which gets less prominent as the T3/V2 cut is tightened. This 
secondary mass peak would be even higher if we loosened the transverse momentum cut 
to pt > 300 GeV, because for moderately boosted top quarks, it is less likely that a single 
jet could capture all three of the subjets coming from the top quark decay. 

In Table [2] and Fig. we see that except for the highest pt range and a too small 
jet radius, a tagging efficiency of around 30% is achievable for a mistagging rate of 1%. 
Note that we do not display results for 300 GeV < px < 400 GeV in Table || nor for 
300 GeV < pt < 450 GeV in Fig. [l^, as in this kinematic range, often not all subjets from 
the decay products of the top quark are captured by the anti-fc^ algorithm, leading to low 
signal efficiencies. This is already noticeable in the 400-500 GeV pt range, where only 45% 
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PT range (GeV) 


400-500 


500-600 


600-700 


700-800 


800-900 


TVo To /To nit, 


.45 


.067 


.58 


.11 


.60 


.13 


.58 


.15 


.53 


.15 


r 3 /r 2 < 0.8 


.43 


.048 


.55 


.080 


.56 


.094 


.54 


.12 


.50 


.12 


t 3 /t 2 < 0.6 


.34 


.016 


.43 


.025 


.44 


.028 


.41 


.039 


.36 


.041 


r 3 /r 2 < 0.4 


.18 


.0019 


.21 


.0028 


.21 


.0035 


.19 


.0038 


.16 


.0069 


1% mistag rate 


.30 


.010 


.34 


.010 


.32 


.010 


.28 


.010 


.21 


.010 


40% tag efficiency 


.40 


.035 


.40 


.018 


.40 


.020 


.40 


.035 


.40 


.061 



Table 2: Tagging efficiencies vs. mistagging rates for top jets : QCD jets with R = 0.8. The top 
row corresponds to just applying the m t0 p invariant mass window (145 GeV to 205 GeV) criteria, 
and the subsequent rows include an additional T3/T2 cut. The bottom two rows indicate the tagging 
efficiencies achievable with a fixed mistagging rate of 1%, and the mistagging rate achievable with 
a fixed tagging efficiency of 40%. 



of the jets pass the mass criterion, whereas the same is true for 58% of the jets in the 500- 
600 GeV pt range. For the same reason, the tagging efficiency is much lower for R = 0.6 



as compared to other jet radii, as seen in Fig. 12(b) . Our method is most effective for top 
jets with 500 GeV < px < 700 GeV. The worse efficiencies for px > 700 GeV are due to 
several factors. A larger portion of top jets have a mass above their mass window due to 
additional FSR at high px, while more light quark and gluon jets will acquire masses high 
enough to pass the mass window cut. Our iV-subjettiness cuts are also less effective at very 
high transverse momenta, since the distinction between 2-subjettiness and 3-subjettiness 



becomes less clear when the three subjets have considerable overlap. Fig. 12(c) compares 
the iV-subjettiness T3/V2 cut to a 2/34/2/23 cut as well as to a benchmark JHTT point, and we 
see that iV-subjettiness compares favorably to these methods for boosted top identification. 
In Fig. we see that a factor of 20 improvement in S/B and a factor of 2 improvement 
in S/y/B is possible using T3/T2 for top jets. 



3.4 Potential Optimization 

In the above efficiency studies, we only considered cuts on the tn/tn-i ratio to tag iV-prong 
boosted objects. One could certainly generalize this approach to include more complicated 
cuts in the tm-\-tn plane or to use a multivariate analysis on the full set of tm values. Such 
studies are beyond the scope of the present work, but as a small step towards optimization, 
we tested the next simplest generalization of tn/tn-\\ a general linear cut in the tn-i~tn 
plane. 

For iV-subjettiness we tested a cut of tn/(tn^\ — rjv-i ; o) < s with the slope parameter 
s <E [0,5]. For YSplitter, we tested the analogous cut on yN,N+i/(UN~i,N — yN-i,N-,o) < s. 
These are essentially the same linear cuts as in the previous subsection except with a 



shift in one of the two iV-subjettiness (or YSplitter) values. In Fig. |9(c)| for W jets, the 



dark grey envelope corresponds to the range n ; o G [—0.075,0.075] with the solid purple 
curve for the special case T2-o = 0. Similarly, the light grey envelope corresponds to 
the range yi2-,o 6 [—25 GeV, 25 GeV] with the dashed purple curve for the special case 
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Top jets vs. QCD jets Top jets vs. QCD jets, 450 GeV < p f < 600 GeV Top jets vs. QCD jets, 450 GeV < < 600 GeV 




Mistagging rate Mistagging rate Mistagging rate 



(a) (b) (c) 

Figure 12: Top jet signal efficiency/background rejection plots. Here, the T3/T2 cut is varied for 
different (a) jet transverse momenta and (b) jet radii. Unless otherwise indicated, the jets have 
R = 0.8 and 450 GeV < pr < 600 GeV. The rightmost points in each plot are for the mt op 
invariant mass window criterium alone, and points to the left of these are obtained with additional 
cuts on the T3/T2 ratio. Figure (c) shows a naive comparison against the YSplitter method, where 
the red lines indicate varying cuts on T3/T2 and 2/34/2/23, respectively, while the shaded bands 
indicate the modified linear cuts of Sec. [TJ. The circle indicates a benchmark JHTT point. In this 
non-optimized test, iV-subjettiness compares favorably to previous top tagging methods. Note the 
different horizontal scale in (c). 



Top jets vs. QCD jets Top jets vs. QCD jets 




(a) (b) 

Figure 13: Improvement compared to an invariant mass cut alone for top jets in (a) signal over 
background (S/B) and (b) signal over square- root background (S/y/B) by using the T3/V2 cut. The 
value of unity on the right corresponds to using only the mt op invariant mass window. 



2/23,0 = 0. The same logic holds for top jets in Fig. 12(c)| with the ranges T2 ; o G [—0.05, 0.05] 
and j/23;0 G [—10 GeV, 10 GeV]. For both methods, small improvements are possible, 
warranting a further exploration of multivariate methods. 



4. Applications to New Physics Searches 

Having seen gains in single W jet and top jet identification in the previous section, we 
now apply iV-subjettiness to the task of heavy resonance reconstruction. Consider a hy- 
pothetical Z' spin-1 resonance which decays to pairs of boosted objects as Z' — > W + W~ 
or Z' — > tt. We will show that an iV-subjettiness cut can lower the minimum detectable 
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cross section for this resonance by up to an order of magnitude with 1 fb _1 of 7 TeV LHC 
data. This illustrates the versatility and power of the iV-subjettiness tagging method in 
the search for new physics. 

4.1 Analysis Overview 



We use the same basic Pythia settings as in Sec. 3.1, simulating pp collisions at y/s = 7 
TeV using the same virtual calorimeter setup. The jet clustering procedure is exactly as 
before, though for all samples we set the anti-/cr jet radius at R = 0.8, which slightly 
degrades the Z' — > W + W~ reach. We perform our selection criteria on events with at least 
two boosted central jets (j>t > 300 GeV and |r/| < 1.3). We impose the same W mass and 
top mass windows and use Tijr\ and T3/T2 as the respective discriminating variables. 

The signal source is a hypothetical "sequential" Z 1 boson of the "extended gauge 
model" [pO|l , one of the reference spin-1 resonance options available in Pythia. We choose 
this type of Z' for convenience since it has couplings both to W bosons and to top quarks. 4 
At the end of the day, we will report the reach in terms of a x Br, so our analysis is 
roughly independent of signal source, apart from parton distribution functions and angular 
correlations, which somewhat affect the signal acceptance. We consider Z' masses between 
750 GeV and 2 TeV, and consider the signal region to be where the combined invariant 
mass of the boosted jets are within 100 GeV of the Z' boson mass. 

There are two main types of backgrounds: reducible backgrounds from QCD dijets and 
irreducible backgrounds from standard model W + W~ and ti production. For simplicity, 
we will only consider QCD dijets for discussing the reach, as the W + W~ and ti processes 
have lower cross sections than the achievable reach with 1 fb" 1 of 7 TeV LHC data. 5 There 
is a potentially important background from VF/Z+jets, where the W/Z decays hadronically 
to form a real boosted W/Z and a jet fakes a boosted W. However, while the W/Z+]ets 
background (with only one jet faking a W) is reduced less by our tagging method than the 
QCD dijet background (with both jets faking Ws), the contribution from VF/Z+jets still 
ends up being about an order of magnitude smaller than that from QCD dijets after our 
optimal cuts. 6 There is also an interesting background from VF-strahhmg that can mimic 
the top jet signal [^4|] which we will not include. 

4.2 Di-W and Di-Top Resonance Results 



To see the effect of the A-subjettiness cut on resonance reconstruction, consider Fig. 14 



for Z — > W + W~ and Fig. 15 for Z' — > ti, where both jets have been tagged as boosted 
hadronic objects. Even for a moderate cut on tjv/tjv-i, there is a substantial decrease in 
the QCD dijet background with only a slight decrease in the resonance signal. 

We can isolate the Z' resonance region by considering dijet invariant masses satisfying 



mz'\ < 100 GeV. In Fig. 16 and Fig. 17, we plot the improvements in S/B and 



4 Of course, for a true sequential Z' , one would likely use the Z' — } £ + £~ mode for discovery. 

5 For ti, we will see that this statement is borderline, so strictly speaking the standard model ti back- 
ground should not be ignored. 

Note that if the mistagging rate is reduced even further with more aggressive iV-subjettiness cuts or 
with an additional selection criteria, it may no longer be possible to ignore the background from W/Z+]ets. 
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Figure 14: Reconstructed invariant mass distributions for (a) a narrow 1 TeV Z' decaying to two 
W jets and (b) dijet QCD events faking Z' — > W + W~ . In both cases, the hardest two jets in the 
event were required to satisfy pr > 300 GeV, \r}\ < 1.3, and 65 GeV < mj e t < 95 GeV, and the 
t~i/t\ criterium is applied to both jets. Note that we can use milder iV-subjettiness cuts than in 



Sec. 3.3, since the cuts are applied to both jets. 
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Figure 15: Reconstructed invariant mass distributions for (a) a narrow 1 TeV Z' decaying to two 
top jets and (b) dijet QCD events faking Z' —> tt. In both cases, the hardest two jets in the event 
were required to satisfy pr > 300 GeV, \rj\ < 1.3, and 145 GeV < mj e t < 205 GeV, and the t 3 /t 2 
criterium is applied to both jets. Again note the milder iV-subjettiness cuts compared to Sec. 



S/yB for &\-W and di-top resonances for three values of the Z' boson mass. Compared to 
the single object efficiencies from Sec. |3], the improvement seen in resonance reconstruction 
approximately factorizes. That is, the tat/ttv-i values for the two hardest jets in any event 
are roughly independent of each other, such that an S/B improvement of e for a given 
Tn/tn-i cut in Sec. [T3] yields an S/B improvement of roughly e 2 for the Z' resonance. 

To see how these improvements in S/B and S/y/B work in practice, we now calculate 
how much an iV-subjettiness cut would lower the minimum cross section for detecting a Z' 
boson in the boosted W + W~ or boosted it channel. For this purpose, we define the fiducial 
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Figure 16: Improvement in (a) S/B and (b) S/y/B for identifying a Z' boson decaying to W 
jet pairs. Compared to the m\y invariant mass cut alone (rightmost points on the plot), cuts on 
./V-subjettiness can significantly improve signal efficiency /background rejection. 
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Figure 17: Improvement in (a) S/B and (b) S/y/B for identifying a Z' boson decaying to top 
jet pairs. Compared to the mt op invariant mass cut alone (rightmost points on the plot), cuts on 
iV-subjettiness can significantly improve signal efficiency /background rejection. 



Z' detection criteria to be twofold. First, we require at least 10 reconstructed candidate 
Z' resonances (S > 10) at a luminosity of 1 fb . Second, we require 5/Err(5) > 5 
("five sigma discovery"), where the combined statistical and systematic error on the QCD 

background is estimated as Err(5) = \J (\/B) 2 + (0.1.B) 2 . 7 Both constraints restrict the 
detectable physical cross section values for the Z' resonance, and we optimize the N- 
subjettiness cuts to satisfy both constraints. 

In Fig. |18(a) , we demonstrate the interplay between the S > 10 and SfErr(B) > 5 
requirements for selecting the optimal T2/T1 cut on the W jets. In Fig. |l8(b)| , we plot the 
lower bound on the detectable cross section for Z', given in terms of cr(Z') x Br(Z' — > 
W + W~) x Br(VF — > jj) 2 . For reference, we also plot the a x Br for the sequential Z' and 
standard model W + W~ production. The T2/T1 cut improves the discovery reach upwards of 
an order of magnitude compared to only using invariant mass to tag W jets. The analogous 
bounds for Z' — > it appear in Fig. 19 where the reach is given in terms of cr(Z') x Br(Z' — > 



7 A 10% systematic on a jet measurement is likely optimistic, but fine as a benchmark for comparison. 
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(a) (b) 

Figure 18: (a): Optimizing the tz/tx cut to determine the minimum physical cross section for 
detection of Z' — > W + W~ with a Z' mass of 1.25 TeV. The value of the cut is chosen such that 
both S > 10 and S/Ett(B) > 5. (b): Lower bounds on the detectable Z' cross section times the 
branching ratio to hadronically-decaying W bosons as a function of the mass of the Z' . Compared 
to a jet invariant mass cut alone, the ./V-subjettiness method to tag boosted Ws gives a substantial 
improvement in the reach. Shown for reference are the a x Br for the sequential Z' model and 
standard model W + W~ production where the di-W invariant mass is within 100 GeV of the 
fiducial Z' mass. 




(a) (b) 

Figure 19: (a): Optimizing the T3/T2 cut to determine the minimum physical cross section for 
detection of Z 1 — > tt with a Z' mass of 1 TeV. The value of the cut is chosen such that both S > 10 
and S/Eti(B) > 5. (b): Lower bounds on the detectable Z' cross section times the branching ratio 
to hadronically-decaying top quark as a function of the mass of the Z'. Compared to a jet invariant 
mass cut alone, the iV-subjettiness method to tag boosted tops gives a substantial improvement in 
the reach. Shown for reference are the a x Br for the sequential Z' model and standard model it 
production where the di-top invariant mass is within 100 GeV of the fiducial Z' mass. 

tt) x Br(t -)• bjj) 2 . The r 3 /r 2 cut can bring Z' signal detection threshold down almost to 
the irreducible ti background level. 8 Thus, iV-subjettiness provides substantial gains for 
discovering new physics with boosted hadronic objects. 

8 The reach for the 750 GeV resonance could likely be improved by using a larger initial jet radius. 
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5. Conclusion 



In this paper, we introduced an inclusive jet shape iV-subjettiness designed to tag boosted 
hadronic objects. We found that the ratio T2/T1 is an effective discriminating variable to 
isolate boosted hadronic W, Z, and Higgs bosons from the background of QCD jets with 
large invariant mass. Similarly, T3/T2 is an effective variable for identifying boosted top 
quarks. As a case study, we observed that A-subjettiness offers significant improvements in 
the detection sensitivity of hypothetical heavy resonances decaying to pairs of electroweak 
bosons or top quarks. Overall, TV-sub jettiness selection methods are at least as good as 
other commonly used discriminating methods for identification of boosted objects. 

./V-subjettiness exhibits several desirable properties which warrant further experimental 
and theoretical investigations. On the experimental side, tn can be calculated on a jet- 
by-jet basis and thereby offers considerable flexibility of application. While we focused 
just on ratios of tn as discriminating variables, multivariate optimization along the lines 
of Sec. 3.4 could improve signal efficiency and background rejection. In addition, some of 
the iV-subjettiness variations mentioned in App. |A| might also be effective discriminating 
variables by themselves or in combination. On the theoretical side, tn is an infrared and 
collinear safe inclusive jet shape which in principle can be defined without the need for 
an algorithmic subjet finding method. Thus, the prospects for theoretical calculations 
involving iV-subjettiness look promising both using fixed-order perturbative calculations 
and using resummation techniques. 

With the first LHC data on the books, the search for new physics is already underway. 
New phenomena may be revealed in the production of highly boosted electroweak bosons 
and top quarks, and we expect that ./V-subjettiness will prove to be a useful variable for 
exploring such extreme kinematic regimes. 



Note Added: While this paper was being completed, Ref. (4^| appeared which defines a 
Lorentz-invariant version of A-subjettiness and uses T2 for boosted Higgs identification. 
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A. Definition of iV-subjettiness 



The definition of iV-subjettiness in Eq. ( |2.lD is not unique, and different choices for rjy 
can be used to give different weights to the emissions within a jet. These generalizations 
of iV-subjettiness are similar to different "angularities" used in e + e~ — > hadrons mea- 
surements. 

Analogously to Ref. p8|, a general iV-subjettiness measure is 



r IT = jy]minMpj,Pt)} , (A.l) 



where d$ is a normalization factor, J runs over the N candidate subjets, and d(pj,pi t ) is 
distance measure between a candidate subjet p j and a jet constituent pk- Like in Sec. \2.2[ 
one needs a method to figure out the candidate subjet directions, which could be achieved 
through a separate subjet finding algorithm or by minimizing tn over possible candidate 
subjets pj. 

There are many choices for d(pj,pk), but a nice two-parameter, boost-invariant choice 
for the distance measure is 

d a 'P(pj,Pk) = PT,k (pt,jT (&Rj,kf • (A.2) 



If desired, one could replace pr,j with £t,j = \Jpt j + m j to include information about 
the subjet mass. 9 For e + e~ applications, one would replace the transverse momentum pr 
with the total momentum \p\ (or the energy E) and AR with the opening angle AQ. A 
natural choice for the normalization factor to keep < r^r < 1 is 

d = max{(p TjJ ) a } (Rof^2pT,k, (A. 3) 

J k 

where i?o is the characteristic jet radius. 

By making d(pj,pk) linear in pT,k, T N is automatically an infrared-safe observable. 
Collinear-safety requires linearity in px,k a s well, but imposes the addition requirement 
that f3 > 0. The value of a is unconstrained. Of course, we are assuming that the 
candidate subjet finding method is also infrared- and collinear-safe. 

In the body of the paper, we used a = 0, j3 = 1. This choice corresponds to treating 
each subjet democratically, and using a fc^-like distance measure. This distance measure 



makes r^v similar to jet broadening [45|, and we found that this was an effective choice for 



boosted object identification. By varying /3, we can change the angular weighting. A thrust- 



like [ 32 1 weighting corresponds to /3 = 2, while other angularities [44] with — oo < a < 2 
are given by f3 = 2 — a. By varying a, we can weight the distance measure by the hardness 
of the subjet directions. Large positive (negative) a means that the minimum in Eq. ( |A.l ) 



9 Obviously, one could also use Ex,k to include the mass of the jet constituent, though in our studies, 
the four- vectors of the calorimeter cells were massless by definition. 

10 By similar, we mean the distance measure has the same Ai?^,fe — ► limit. Because thrust-like ob- 
servables are defined in a preferred rest frame and we are working with a longitudinally boost-invariant 
measure, the correspondence is inexact. 
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Figure 20: Top row: W jets. Bottom row: QCD jets with invariant mass close to mw The 
coloring and labeling is the same as in Fig. [p. The title of the plot gives the calculated value of 
T2/TI7 and generically W jets have smaller T2/T1 values than QCD jets of comparable invariant 
mass. 



is given by the distance to hardest (softest) candidate subjet. Further studies of boosted 
object identification using different values of a and (5 would be interesting, since studies 
of jet angularities have shown that additional information about jet substructure can be 
gleaned by combining different angular information [29]. 

B. Additional Event Displays 

To further demonstrate our subjet finding method and the distinguishing power of N- 



subjettiness, we compare W jets and QCD jets in Fig. |20| and compare top jets and QCD 



jets in Fig. 21. The discriminating variable tjv/tjv_i measures the degree to which the jet 
energy is aligned along the ./V candidate subjet directions compared to the N — 1 candidate 
subjet directions. 
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Figure 21: Top row: top jets. Bottom row: QCD jets with invariant mass close to m top - The 
coloring and labeling is the same as in Fig. |[ The title of the plot gives the calculated value of 
T zl T 2i an d generically top jets have smaller t 3 /t 2 values than QCD jets of comparable invariant 
mass. 
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Figure 22: Distributions of (a) 2/12 and (b) 2/23 for boosted W and QCD jets. For these plots, we 
impose an invariant mass window of 65 GeV < mj e t < 95 GeV on jets of R = 0.6, px > 300 GeV, 
and \r)\ < 1.3. 
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Figure 23: (a): Distribution of 1/23/2/12 for boosted W and QCD jets, (b): Density plot in the 
2/12-2/23 plane. 



C. Comparison to YSplitter 



In Sec. IH, we compared iV-subjettiness to a naive application of the YSplitter method 
H ||, U|. Here, we collect various plots of yN,N+i and their ratios so the reader can visually 
compare the discriminating power of iV-subjettiness and YSplitter. The results for W jets 

As mentioned in the text, 



are shown in Figs. |22| and|23|, and top jets in Figs. |25, and|2f^ 
a full comparison of the two methods would require an optimization of all cuts, though it is 
encouraging that in Sec. we found that linear cuts in the t^-t^^i plane are generically 
more effective than linear cuts in the yN,N+\-VN-i,N plane. 
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Figure 24: Distributions of (a) 2/12, (b) 2/23, and (c) 2/34 for boosted top and QCD jets. For these 
plots, we impose an invariant mass window of 145 GeV < m- ict < 205 GeV on jets with R = 0.8, 
p T > 300 GeV, and \t}\ < 1.3. 



145 GeV<m < 205 GeV 



0.06r 

0.05 

CD 

| 0.04 

o 0.03 
<u 

| 0.02 
0.01 












Top jets 

QCD jets 













































0.2 0.4 0.6 0.E 

y 23 /y 12 ofjet 
(a) 



145 GeV<m < 205 GeV 



0.08 



o 0.04 



^ 0.02 









Top jets 

QCD jets 







































0.2 0.4 0..6 0.£ 



(b) 



Figure 25: Distributions of (a) 2/23/2/12 and (b) 2/34/2/23 for boosted top and QCD jets. Analogous 
to Fig. t^, we do not use 2/23/2/12 for top tagging in this paper, though it does contain additional 
information. 
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Figure 26: Density plots in the (a) 2/12-2/23 plane and (b) 2/23~2/34 plane for boosted top and QCD 
jets. 
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